Inverse problem: Reconstruction of modified gravity action in Palatini formalism by 

Supernova Type la data 
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We introduce in f(R) gravity-Palatini formalism the method of inverse problem to extract the 
action from the expansion history of the universe. First, we use an ansatz for the scale factor and 
apply the inverse method to derive an appropriate action for the gravity. In the second step we 
use the Supernova Type la data set from the Union sample and obtain a smoothed function for the 
Hubble parameter up to the redshift 1.7. We apply the smoothed Hubble parameter in the inverse 
approach and reconstruct the corresponding action in f(R) gravity. In the next step we investigate 
the viability of reconstruction method, doing a Monte-Carlo simulation we generate synthetic SNIa 
data with the quality of union sample and show that roughly more than 1500 SNIa data is essential 
to reconstruct correct action. Finally with the enough SNIa data, we propose two diagnosis in 
order to distinguish between the ACDM model and an alternative theory for the acceleration of the 
universe. 

PACS numbers: 04.50. +h, 95.36.-f-x, 98.80.-k 
I. INTRODUCTION 



Combination of SNIa+CMB data shows that universe 
is in the positive acceleration phase [1,2]. This result is 
in contradiction with our expectation from the behavior 
of the ordinary matter. The simplest solution is assum- 
ing a cosmological constant in the Einstein field equation 
as a constant of integration [3] . While this simple modi- 
fication explains the observational data [4], however the 
cosmological constant suffers from the fine tuning and 
coincidence problems [5]. One of the solutions is intro- 
ducing a scalar field which provides a time dependent 
negative equation of state [6] . 

The other possibility is the modification of the grav- 
ity law in such a way that it behaves as standard Gen- 
eral Relativity in strong gravitational regimes and repulse 
particles in the low density cosmological scales [7]. The 
modified gravity models can be examined with three cat- 
egory of observations of (a) cosmological dynamics [8], 
(b) local gravity [9] and (c) the evolution of large scale 
structure [10]. 

In this work we use SNIa data to reconstruct an appro- 
priate f(R) gravity model in Palatini formalism [11] with 
the inverse method. The method is the extension of the 
work introduced for the metric formalism by Rahvar and 
Sobouti in [12]. The inverse method also is introduced 
in the work by Capozziello et al. [13] in the metric for- 
malism. In this method we need to know the dynamics 
of Hubble parameter from the observational data. Many 
methods for the extraction of the Hubble parameter have 
been introduced in the literature [14]. Here we use the 
smoothing method suggested by Shafieloo et al. [15] and 
apply it to the SNIa Union sample [16]. We use H(z) in 
the inverse method algorithm to reconstruct the corre- 
sponding action. We test the reliability of this method by 
doing a Monte-Carlo simulation and generating the syn- 



thetic SNIa data according to an action and comparing 
the reconstructed action with the original one. Finally 
we introduces two diagnosis for distinguishing standard 
ACMD model from the alternative models. 

The structure of this article is as follows: In section 
II we introduce f{R) modified gravity in Palatini for- 
malism, derive the equation of motion and obtain the 
dependence of Hubble parameter to the Ricci scalar and 
scale factor. In section III we introduce the method of 
inverse problem in Palatini f(R) gravity. In section IV 
we use the method to the real data of SNIa, smoothing 
the supernova data we extract the Hubble parameter in 
terms of redshift and apply it to extract the action [15]. 
Also in order to show the level of confidence of our re- 
sults, we simulate 100 realization of SNIa data to extract 
the Hubble parameter and compare it with that obtained 
directly from fitting to the model. In the section V, we 
examine the viability of reconstruction method and de- 
pendence of the results to the number of SNIa data. In 
section VI we propose two diagnosis as probe to distin- 
guish between the ACDM and the alternative models. 
Section VII concludes the paper. 



II. MODIFIED GRAVITY IN PALATINI 
FORMALISM 

For f(R) gravity, there are two main approaches to ob- 
tain the field equation. The first one is so-called metric 
formalism which is obtained by the variation of the action 
with respect to the metric. In this case the derived field 
equation is a fourth order nonlinear differential equation. 
In the second approach, which is called Palatini formal- 
ism, the connection and metric are considered as inde- 
pendent fields and the variation of action with respect 
to these fields results in a set of second order differential 
equations. The Palatini formalism is a plausible candi- 



date to be the effective classicai theory of gravity from a 
more fundamental theory of Loop Quantum Gravity [17]. 

Let us take a general form of the action in Palatini 
formalism as: 



S[f;g,T,* r . 



h 1 ^ 



(i) 



where n = 8itG and S m [g M „ , $ ra ] is the action of mat- 
ter depends on the metric g^v and the matter field 
R = R(g,T) = g^R^ip) is the generalized Ricci scalar 
and R^ v is the Ricci tensor made of affine connection. 
Varying action with respect to the metric results in: 



f(R)R flv (t) - l -f{R)g^ = kT^, 



(2) 



where prime is the derivative with respect to the Ricci 
scalar and T^„ is the energy momentum tensor 
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(3) 



On the other hand varying the action with respect to the 
connection results in: 



V a [f(R)V^g 



inn 



0. 



(4) 



where V is the covariant derivative defined from par- 
allel transformation and is given by the affine con- 
nection. From equation (4), we define a new metric, 
= f'{R)g^y conformally related to the physical met- 
ric where the connection is the Christoffel symbol of this 
new metric. 

We apply flat FRW metric (namely K — 0) for the 



universe 



ds 2 = -dt 2 + a(t) 2 8 i jdx i dx j , 



(5) 



and assume that universe is filled with a perfect fluid with 
the energy-momentum tensor of = diag(—p,p,p,p). 
Using the metric and energy momentum tensor in Eq.(2), 
we obtain the generalized FRW equations. It should be 
noted that the conservation law of energy-momentum 
tensor, T^ v ;ll = is satisfied according to the covari- 
ant derivative with respect to the metric and this defini- 
tion guarantees the motion of particles on geodesic [18]. 
Combination of Goo and G\ results in [19]: 



(H 
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Taking the trace of equation (2) results in 
Rf'(R) - 2f(R) = kT, 



(6) 



(7) 



where T = g^T^ = —p + 3p. The time derivative of this 
equation results R in terms of time derivative of density 



and pressure. Using the equation of state of cosmic fluid 
p = p(p) and the continuity equation, the time derivative 
of Ricci obtain as: 



R = 3kH 



(l-3dp/d P )(p + p) 
Rf" - f'(R) 



(8) 



To obtain generalized first FRW equation we start with 
equation (7) and obtain the density of matter in terms 
of Ricci scalar as: 



up 
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where w = p/p. We substitute equation (9) in (6) and 
use equation (8) to change the time derivative to d/dR. 
We rewrite equation (6) as follows: 



H 2 = 



1 
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(10) 



On the other hand using equation (7) and continuity 
equation, the scale factor can be obtained in terms of 
Ricci scalar 



1 



npo(l — 3uj) 



(2/ - Rf) 



(11) 



where po is the energy density and ao is the scale factor 
(set to one, i.e. ao = 1) at the present time. Now for a 
generic modified action, omitting Ricci scalar in favor of 
the scale factor between equations (10) and (11) we can 
obtain the dynamics of universe (i.e. H = H(a)). 

For the simple case of matter dominant epoch uj = 
which is in our concern, these equations reduce to: 
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INVERSE METHOD IN PALATINI 
FORMALISM 
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(13) 



In this section we introduce the inverse method to ex- 
tract f(R) action in Palatini formalism from the dynam- 
ics of the universe. This method has been studied in the 
work by Rahvar and Sobuti in the metric formalism [12] 
and we extend it to the Palatini formalism. 

Replacing / with the first derivatives of action from 
the equation (7), 



f(R) = ±[RF-K(3p-p)}, 



(14) 



equation (6) can be written as follows: 



1F N , R k , 



(15) 



where F is defined as F = df/dR. It should be noted that 
Ricci tensor in the Palatini formalism is given in terms 
of the conformal metric of h^ v = F(R)g^. Substituting 
the new metric in the definition of Ricci scalar results in: 



R^v = R^{g) 



3V M FV„F 
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where R fiu (g) is the Ricci tensor defined in terms of met- 
ric gp,v By taking trace from equation (16) we obtain 
relation between the Ricci scalar in Palatini and metric 
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On the other hand using the FRW metric, R(g) as the 
Ricci scalar in metric formalism is given by the Hubble 
parameter as: 



R{g) = 6H+12H 2 



(18) 



where for simplicity in the calculation we rewrite this 
equation by changing the time derivative to the rcdshift 
derivative. In what follows we use prime for the deriva- 
tive with respect to the redshift. 



R(g) = -6HH'(1 + z) + 12H 2 



(19) 



For the Ricci scalar in Palatini formalism from equa- 
tion (17), we obtain: 



R = -6HH'(1 + 
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FIG. 1. Deviation parameter £ = F(R) — 1 is obtained 
from the ansatz scale factor (solid-line) is compared with 
^-parameter of 1/R modified gravity (dashed-line) in terms 
of redshift. For the redshifts z > 1, the f{R) gravity model 
converge to the Einstein-Hilbert action. 



On the other side we can calculate the Ricci scalar from 
equation (20) in terms of redshift. Eliminate z in favor 
of Ricci scalar results in F(R). Finally by the numerical 
integration of this function we can obtain the modified 
gravity action, f(R). 

In the rest of this section to test this algorithm we 
use an ansatz for the scale factor and try to extract the 
corresponding action from a given dynamics. We apply 
the following ansatz for the scale factor proposed in [12]: 



a(r) 



1 



1 
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(22) 



Substituting this equation in (15) results in a differ- 
ential equation for the evolution of F as a function of 
redshift: 



3F'' 



F - -- 
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H(l + z) H 2 {l + z)' 



;(p + p)=Q- 



(21 



In the matter dominant epoch we rewrite this equation 
by putting p = and p = po(l + z) 3 , where we can replace 
upo with 3-ffgf2 m . It should be noted that the definition 
of the O m is different from that in the standard FRW 
equations. 

For a given dynamics of universe (i.e. a(t)), we can ex- 
tract the Hubble parameter in terms of redshift, applying 
it in equation (21) will provide F in terms of redshift. 



(^'in which t = tH$ is a dimensionless time parameter de- 
fined in the interval of to 1 and Hq is the Hubble param- 
eter at the present time. This proposed dynamics has two 
free parameters of a and p. We obtain the corresponding 
Hubble parameter from this scale factor and consequently 
the distance modulus and compare the model with the 
observed SNIa data. The best value for the parameters 
of model has been obtained roughly, p = ^ and a = 2 
in [12]. Substituting the scale factor in equation (21), 
we obtain the deviation parameter of £ = F — 1 from 
GR as a function of redshift, which is plotted in Figure 
(1). In order to reconstruct the action in terms of Ricci 
action, we obtain the Ricci scalar in terms of rcdshift 
and finally by eliminating rcdshift between F and Ricci 
scalar, we obtain F = F(R). Integrating this function 
provides action in terms of Ricci scalar. We also plot 
the ^-parameter of 1/R modified gravity for comparison 
with the reconstructed modified gravity in Figure (1) and 



0.22 



0.2 



Reconstructed from Ansatz 
1 R Modified gravity 




FIG. 2. The difference of modified gravity action from the 
Einstein-Hilbert action is plotted versus Ricci scalar (solid- 

4 

line) and compared with the action of f(R) = R— *g- (dashed- 
line) . The Ricci scalar is normalized to its value at the present 
time, R„. 



R — f(R) as a function of R in Figure (2). Comparison of 
the extra term with the Einstein-Hilbert action roughly 
resembles to /U 4 /i? function with /i 2 ~ 10 _1 i? p where R p 
is the present value of Ricci scalar. 



IV. RECONSTRUCTION OF THE DYNAMIC OF 
THE UNIVERSE BY SNIA DATA 

In section III, wc showed that it is possible to recon- 
struct the modified gravity action by knowing the dy- 
namics of the universe. In this section, we use SNIa 
cosmological data to obtain the Hubble parameter and 
consequently reconstruct the action of modified gravity. 

The dynamics of Hubble parameter, H(z), can be ob- 
tained if the distance modulus of Supernovas data as a 
function of redshift is known. In FRW-fiat universe, the 
Hubble parameter is related to the distance modulus of 
SNIa as follows: 



(23) 



The main challenging point in this procedure is the lim- 
ited number of observed SNIa which impose an uncer- 
tainty in calculating the continues function for H(z). The 
overall number of Supernovas which has been detected is 
in the order of 300-400. We use the latest Supernova 
data of the Union sample to extract the Hubble param- 
eter [16]. To make a continues Hubble parameter, we 
follow the procedure known as reconstruction method, 
proposed by Shafieloo et al. in [15]. In this algorithm, 



a non-parametric function is used for smoothing the dis- 
tance modulus of supernova data over the redshift. Here 
we choose a guess model resemble to the observed dis- 
tance modulus of the supernova data. In the next step we 
subtract the distance modulus of the observed data from 
the guess model using a gaussian function for smoothing 
the observed data as follows: 



i obs ( Zl )} exp- 



2A 2 



where 



iV(z)- 1 =^exp[ 



2A 2 



(24) 



(25) 



where Zi represents the redshift of each SNIa in the Union 
sample. The sum term is considered for all 307 Union 
sample data, fi obs (zi) is the observed distance modulus 
and fi 9 (z) is a continues guess model for the distance 
modulus. N(z) is a normalization factor and A is a suit- 
able redshift window function. fi d (z) is a smoothed con- 
tinues function for the residual of luminosity distance in 
terms of redshift. Now the corrected distance modulus 
is added to the guess model to generate the new smooth 
function for the distance modulus: 



^(z) = ^(z) +t i d (z). 



(26) 



We repeat this procedure using [i s {z) as the new guess 
function. It can be shown that after a finite time of this 
irritation, \ 2 of the smoothed function with respect to 
observed data will converge to a fixed value. It means 
that we will find a continues distance modulus function 
with the best fit to the real data. It is shown in [15] 
that the result of the best continues function is indepen- 
dent of the choice of the first guess model. Having the 
the smoothed luminosity distance we use equation (23) 
to obtain the Hubble parameter H(z). Another point is 
that the results will clearly depend upon the value of A 
in equation (24). A large value of A produce a smooth 
result, but the accuracy of reconstruction worsens, while 
small value A gives more accurate, but noisy result. Con- 
sidering the frequency of SNIa data observed in the Union 
sample, we choose A = \/a7(1 + z)Aq for the window 
function, where M is the number of irritation need to 
converge x 2 and Ao = TV -1 / 3 where N is the number of 
SNIa [15]. 

In what follows we find the uncertainty in H (z) from 
this method to use it for reconstructing the appropriate 
action in the Palatini formalism. For this purpose we do 
a Monte-Carlo simulation, generating 100 realization of 
SNIa data and using the same distribution of SNIa in 
terms of redshift reported by Kowalski et al. [16]. Also 
we use the same error bars of distance modulus in the ob- 
served data. In order to simulate the synthetic distance 
modulus of SNIa, wc assume a dark energy model for the 
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FIG. 3. The uncertainty in E(z) = -jj- from the smoothing 
method, resulting from the Monte-Carlo simulation. We gen- 
erate one hundred realization for the distance modulus and 
obtain the Hubble parameter which all these continues func- 
tions reside inside the boundaries in the figure. 



universe with a constant equation of state of lu = —0.75. 
Choosing ACDM as the guess model for these data, we 
obtain the Hubble parameter E{z) = H{z)/Hq for 100 re- 
alization of supernova data. Figure (3) shows the bound- 
aries for E(z), resulted from 100 realization of Supernova 
data. 

Now we can apply this uncertainty to H(z), resulting 
from the smoothing procedure on the real Union sam- 
ple [20]. We plot the Hubble parameter from the Union 
sample as shown in Figure (4) with a margin represents 
the uncertainty that is obtained form the Monte-Carlo 
simulation. Using the method described in section III, 
we extract the corresponding modified gravity in Pala- 
tini formalism. The parameter of £ = df(R)/dR — 1 in 
terms of Ricci scaler shows a small deviation from the 
Einstcin-Hilbert action as shown in Figure (5) with the 
uncertainty of this parameter. Here the deviation from 
the Einstein-Hilbert action is in the order of £ ~ 10~ 3 . 
We ask the reliability of this result in terms of the num- 
ber of SNIa data. In the next section we will discus about 
this issue. 

V. VIABILITY OF SMOOTHING METHOD AND 
THE NUMBER OF SNIA DATA 

In this section we examine the viability of smoothing 
of the Hubble parameter. The results will be applied to 
the modified gravity models in the next section. 

According to the methodology of the smoothing pro- 
cedure, the Hubble parameter depends on the quantity 



FIG. 4. Dashed line is the reconstructed Hubble parameter 
from SNIa Union sample data. Bold lines are the confidence 
level of Hubble parameter. 



and quality of SNIa data. Similar to the simulation in 
the previous section we generate 100 realization of 307 
SNIa data, using the redshift and the uncertainty of the 
distance modulus in the Union sample within the frame- 
work of ACDM model. We generate continues Hubble 
parameter with the margin of uncertainty which is shown 
in Figure (6). On the other hand we obtain the Hubble 
parameter in ACDM model though fitting the observed 
distance modulus of the data to the model as shown in 
Figure (6). The margin represents la level of confidence 
for the Hubble parameter. Comparing the uncertainties 
from these two methods show that the smoothing method 
is not reliable with 307 number of supernova data. 

We increase the number of SNIa data with the qual- 
ity of Union sample to have the same uncertainty in the 
Hubble parameter from the two methods. For N > 1500 
as shown in Figure (6) we can rely on the smoothing 
method. We can repeat this simulation with precise data 
than the Union sample to achieve this goal with a smaller 
number of Supernova data. 

VI. DISTINGUISHING BETWEEN THE 
MODIFIED GRAVITY MODELS AND ACDM 

In this section we present two diagnosis in order to dis- 
tinguish between ACDM from an alternative model. The 
first method is (a) Comparison of the Hubble parameters 
and the second method is (b) the S-function. Both of di- 
agnosis are applicable by using the inverse method. 
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FIG. 5. The bold line shows the reconstructed modified 
gravity £ = F — 1 parameter versus Ricci scalar normalized 
to its present value. The dashed lines show the confidence 
level of the £ parameter corresponds to the uncertainty of the 
Hubble parameter resulting from the smoothing procedure. 



FIG. 6. Dashed lines and dotted dashed line are the recon- 
structed Hubble parameters from 307 and 1535 SNIa data, 
respectively. Solid lines are 1-a confidence level of Hubble 
parameter from data fitting to 307 and 1535 data in ACDM 
model. 



A. Comparison of the Hubble parameters 

In the previous section we have seen the smoothing 
method to generate a continues Hubble parameter from 
the observed data. If the dynamics of universe follows 
rather than the ACDM model, can we distinguish the real 
model of the universe, knowing the Hubble parameter 
from the data ? 

We can compare the reconstructed Hubble parameter 
directly from the observed data with that of obtained 
from fitting data to the ACDM model. Wc subtract the 
two Hubble parameters and compare it with la deviation 
of the smoothed Hubble parameter. As we simulated 
in the previous section using N > 1500 is sufficient to 
reliable on the uncertainty of the Hubble parameter. 

4 

We use f{R) = R — ^ action in Palatini formalism 
and generate 1535 supernova data with the quality of 
the Union sample. Using the smoothing method we ob- 
tain the Hubble parameter. On the other hand we fit the 
generated data with the ACDM model. Figure (7) com- 
pares the two Hubble parameters and shows that for the 
redshifts with z > 1, the two models arc distinguishable. 



B. The method of H-function 

We define a new criterion to distinguish the ACDM 
from an alternative model as follows: 



„ 2E(z)E'(z) 

"~ 3(1 + z) 2 ' 10 

where prime is differential with respect to the redshift. 
For the case of F = 1 in ACDM universe, from equation 
(21), S reduces to S = f2^. If we have the smoothed 
normalized Hubble parameter from the observed data, 
E(z), then 2 can be obtained. Any deviation from a con- 
stant value for this function is a diagnosis for the ACDM 
universe. This method is applicable both for the dark 
energy and modified gravity models. Let us take oj(z) as 
the equation of state of dark energy, we can use FRW 
equations to obtain S as follows: 

z = n° m + (i-nl)±±^l, (28) 

(1 + z) A 

where for the constant value of the equation of state, 
ui(z) = — 1, this equation reduces to H = fl^. 

In order to show the deviation from the Einstein- 
Hilbert action, we can calculate S — flm as a function 
of redshift. Here we know from a direct cosmolog- 
ical observations such as gravitational weak lensing. To 
show how this method works we generate SNIa data from 
f(R) = R — fJ. 4 /R action resemble to the Union sample 
and after smoothing the Hubble parameter, we plot S 
as a function of redshift. We generate S function for two 
cases with 307 and 1535 supernova data as represented in 
Figure (8). To compare this function with a direct mea- 
surement of Qm j we use the weak lensing data which 
puts limit on the matter content of the universe in the 




FIG. 7. The solid-line is the Hubble parameter with 1-a er- 
ror bar from fitting the simulated data to the ACDM model. 
The dashed-line is the confidence level of Hubble parameter 
extracted from the smoothing method of the Hubble parame- 
ter. The figure shows the deviation of the two curve for z > 1. 



FIG. 8. S-function versus redshift is plotted for two series 
of SNIa data, with la level of confidence. The dash-lines 
and dashed-dotted lines indicate for 307 and 1535 data re- 
spectively. Solid lines indicate the matter density of universe 
derived from gravitational weak lensing. 



range of Ctm = 0.2±q;| [21]. As shown in Figure (8), in- 
creasing of SNIa data will pin down the 3 function more 
precisely where today's SNIa data is not sufficient. On 
the other hand direct observations of matter density of 
universe are not much precise to claim any model distin- 
guishing results. 



of Supernova data from the other hand will enable us to 
identify either our universe follows ACDM universe or 
some modification is necessary. 



VII. CONCLUSION 

One of the most puzzling questions in the cosmology 
is the physical mechanism for the acceleration of the uni- 
verse. Is it driven by a cosmological constant or universe 
is filled with an exotic dark energy or the Einstein grav- 
ity should be modified? In this work we proposed an 
inverse method to extract the action of a modified grav- 
ity in the Palatini formalism from the expansion history 
of universe. 

We used the smoothing method to obtain a continues 
Hubble parameter from the Supernova Type la Union 
sample data. We showed that more than 1500 supernova 
data with the quality of union sample is essential to have 
a reliable Hubble parameter. Finally we proposed two 
cosmological diagnosis in order to distinguish between 
ACDM and alternative models. The first one compares 
the smoothed Hubble parameter from the SNIa data with 
the standard ACDM model. In the second approach we 
define a new parameter where it is constant, equal to the 
f2 m for the ACDM universe and vary with the redshift for 
any alternative model. A precise measurement of matter 
content of the universe from one hand an enough number 
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